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> ■ ABSTRACT 

a^ 

^ I We consider the competition of the different physical processes that can affect 

the evolution of a flame bubble in a Type la supernovae — burning, turbulence 
and buoyancy. Even in the vigorously turbulent conditions of a convecting white 
dwarf, thermonuclear burning that begins at a point near the center (within 



O 

o 

Qh' 100 km) of the star is dominated by the spherical laminar expansion of the flame, 

O ' until the burning region reaches kilometers in size. Consequently flames that ig- 

^ • nite in the inner ^ 20km promptly burn through the center, and flame bubbles 

■ anywhere must grow quite large — indeed, resolvable by large-scale simulations 
of the global system — for significant motion or deformation occur. As a result, 

^ ' any hot-spot that successfully ignites into a flame can burn a significant amount 

■ of white dwarf material. This potentially increases the stochastic nature of the 
explosion compared to a scenario where a simmering progenitor can have small 
early hot-spots float harmlessly away. Further, the size where the laminar flame 
speed dominates other relevant velocities sets a characteristic scale for fragmen- 
tation of larger flame structures, as nothing — by definition — can easily break the 
burning region into smaller volumes. This makes possible the development of 
semi-analytic descriptions of the earliest phase of the propagation of burning in 
a Type la supernovae, which we present here. Our analysis is supported by fully 
resolved numerical simulations of flame bubbles. 



Subject headings: supernovae: general — white dwarfs - hydrodynamics — nu- 
clear reactions, nucleosynthesis, abundances — conduction — methods: numer- 
ical 



- 2 - 



1. INTRODUCTION 

Numerical modeling of Type la supernovae (SNe la) has shown that a thermonuclear 
burning front propagating outward from the center of a carbon/oxygen white dwarf can 
release enough energy to unbind the star (see for example Ropke and Hillebrandt 2005; 
Gamezo et al. 2005). To prevent the overproduction of iron-group elements and allow for 
the necessary expansion of the white dwarf, this burning front must begin as a deflagration 
(Nomoto et al. 1976), but may later transition into a detonation (Niemeyer and Woosley 
1997; Khokhlov et al. 1997). One of the greatest uncertainties in modeling SNe la lies in 
understanding the earliest stages of burning. Large-scale simulations seed the explosion with 
one or more physically large hot-spots. Variations in the number and location of these hot- 
spots can lead to a wide variety of different outcomes (Niemeyer et al. 1996; Plewa et al. 
2004; Garcfa-Senz and Bravo 2005; Roepke et al. 2006). 

Leading up to the explosion, accretion onto the surface of the white dwarf compresses 
the star and raises its central density and temperature. As the temperature increases, the 
carbon near the center can begin to fuse, releasing heat and driving convection throughout 
the star. This smoldering phase lasts for about a century (Wunsch and Woosley 2004). 
Numerical simulations have just begun to capture the dynamics of this phase (Hoflich and 
Stein 2002; Kuhlen et al. 2006). Some models have shown that the central heating leads to 
a large scale dipole flow in the star (Kuhlen et al. 2006) which may favor off-center ignition. 
What we know for sure is that at the point of ignition, the white dwarf will be turbulent, 
and likely with typical intensities of 100 km s~^ on scales of 200 km (Woosley et al. 2004) 

As the nuclear energy generation rate increases, eventually, hot plumes will develop 
whose cooling time (via expansion and neutrinos) is longer than the nuclear timescale, and a 
flame will be born (Woosley et al. 2004). Once ignited, the flame bubble will burn outward 
as it buoyantly rises, eventually wrinkling due to drag, the Rayleigh- Taylor instability, and 
interactions with turbulence generated from the convection. In this paper, we investigate the 
evolution of these very first flame 'bubbles'. We look at the interplay between these physical 
processes to determine what size the bubble can grow to while remaining roughly spherical — 
this is the maximum size bubble that large scale calculations should use to initialize their 
explosion. We also are interested in the backward propagation of the flame and determining 
how off center a flame can begin such that buoyancy cannot carry it away faster than it 
burns through the center of the star. This will have important consequences for off-center 
ignition models. 

Garcia-Senz and Woosley (1995) explored the evolution of buoyant plumes in the pre- 
ignition white dwarf and argued that multi-point, rather than single point ignition is more 
likely. Their results showed that the hot spots can travel great distances from the center 
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before igniting, and that this distance is relatively insensitive to the initial bubble size. By- 
chkov and Liberman (1995) described a picture of SNe la explosion in which the flame begins 
at multiple points near the center which quickly burn together. The resulting merged bub- 
ble has a higher terminal velocity and floats away from the center, leaving behind unburned 
carbon and oxygen at the center. They argue that this sets the stage for a detonation wave, 
ignited near the center, to burn through the star. 

More recently, lapichino et al. (2006) present a parameter study of two-dimensional 
smoldering bubbles in pre-ignition conditions in the white dwarf. They consider initial sizes 
in the range of 0.2 to 5 km at a range of initial distances from the center of the star. They 
conclude that larger bubbles are more likely to run away and that, in contrast to Garcia- Senz 
and Woosley (1995), the short timescale before these bubbles either ignite or disperse implies 
that the bubbles do not stray far from their origin. They attribute the disparity with the 
Garcia-Senz and Woosley (1995) results to multidimensional effects in the bubble evolution. 

We concern ourselves with a slightly different problem — once a flame ignites at a point 
and begins burning outwards (a flame bubble), how does it evolve. This has seen very limited 
numerical work, mostly because of the difficultly in resolving the flame and evolving the state 
for long periods of time. Initial results of resolved bubbles in three- dimensions were presented 
in Zingale et al. (2005a). Here, we restrict our simulations to two dimensions, to allow for a 
range of different calculations to be carried out while resolving both the flame structure and 
the large-scale hydrodynamic motions. Since we consider only the onset of bubble distortion 
through hydrodynamic motions, two dimensional simulations are quite adequate. 

In § 2, we present a semi-analytic model for the evolution of a flame bubble. Numerical 
simulations supporting this model are presenting in § 3. Finally, in § 4, we discuss the 
implications of this flame bubble model on the ignition of Type la supernovae. 

2. THE EVOLUTION OF A FLAME BUBBLE 

We consider the evolution of a flame bubble by comparing the timescales of the physical 
processes that can affect the dynamics. For SNe la, these timescales are set by the laminar 
flame speed. Si, the bubble rise velocity, Vj,, and the turbulent intensity, Wturb; over scales 
of the size of the bubble and smaller. If the laminar flame speed exceeds all of the other 
velocities on the bubble scale, then we expect the bubble to remain mostly spherical in shape 
and expand in radius solely by burning outward. 

For a spherical flame burning outward — or indeed burning in any geometry where the 
burned material is not free to move — the expansion of the material as it is burned causes the 
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bubble to expand at a rate greater than the laminar flame speed. The mass of fuel consumed 
in a time interval At is AnR'^pf^ciSiAt, where R is the radius of the bubble. Across the burning 
front, the ash expands, pushing the radius of the bubble out at a speed R. The ash mass 
formed over the time interval At is equal to the mass consumed, so 

47ri?Vfuci5iAt = 47ri?Vashi?At (1) 

Therefore, the expansion velocity of the bubble is 

R='-^S,= '-^S, (2) 

Pash ' 1 - At ' ^ ' 

where At = (pfuci — Pash) / (Pfuci + Pash) IS the Atwood number. The radius of the bubble would 
be given by 

1 -I- At 

R{t) = i?o + y^aF^'^ • 

Besides this 'volume creation' effect, geometry also comes in through the curvature 
and strain of the outgoing spherical flame. The curvature and strain of the flame 'dilute' 
heat transport outwards, decreasing the velocity of an outgoing spherical flame {e.g., Dursi 
et al. 2003). This effect is only significant for bubbles sufficiently small (with respect to the 
flame thickness) that the curvature and strain are relatively large, typically tens of flame 
thicknesses. We will generally be considering flame bubbles of at least many thousands of 
flame thicknesses in radius, where this correction is completely unimportant; however we 
include it in the calculations here for consistency when we later consider very low density 
flames where the thickness is large enough that this effect is significant. The correction to 
the flame velocity is 

S'i = Si{l + MaKa) (4) 

where the Markstein number Ma for astrophysical flames is of order —1 (Dursi et al. 2003), 
and the dimensionless strain rate, given by the Karlovitz number Ka, is given for an expand- 
ing flame as 

Ka ^ (Izilill ,5) 

where {d—l)=2 for an expanding spherical flame, and 1 for an expanding flame 'cylinder'. 
Including both the volume creation and curvature effects gives 



R = S, 



1 + At 
1 - At 



, /1 + At 



1 - Ma^{d- 1 



R' ' VI - At 



-1 



(6) 



Competing against the simple expansion of the flame are the buoyancy caused by the 
density difference between the cold fuel and hot burned ash, and externally imposed turbu- 
lence from the centuries-long convective simmering which has occurred before the ignition 
of the first flames. 



- 5 - 



The buoyancy of the bubble causes the bubble both to rise, and to potentially be 
disrupted by the vortical motions generated by its own rising or Rayleigh- Taylor /Kelvin- 
Helmholtz instabilities {e.g., Robinson et al. 2004; lapichino et al. 2006). Buoyancy causes 
the burned material to accelerate upwards from the center of the white dwarf, however very 
quickly the bubble attains a (quasi-) steady terminal velocity set by the balance between 
gravitational acceleration and fluid drag. 

In principle, obtaining this steady characteristic buoyant velocity is almost intractable, 
as the bubble will distort (in a velocity-dependent way) as it rises, greatly complicating 
the flow around the bubble. However, since we are interested in the case where the flame 
velocity is still nearly dominant, we expect the flame bubble to remain essentially spherical. 
In this limit there is extensive theoretical work {e.g., Harper 1972; Magnaudet and Eames 
2000) both past and ongoing {e.g., Joseph 2003). Such work generally considers terrestrial 
fluids, where the gravitational acceleration is spatially constant, and fluid properties do not 
change over the scale of the bubble. For the analysis considered here, we will make the same 
assumptions; in a white dwarf, the typical pressure scale height is something like 200 km — as 
we will show later, we are mostly concerned with bubbles that are < 1 km, so changes in 
gravity or fluid state variables across the size of the bubble are not expected to be important. 

In the case of a bubble rising in a quiescent medium, the drag is due to the wake behind 
the bubble as it rises. A gas bubble {e.g., one where the bubble density is negligible) in 
an incompressible, inviscid fluid, has been shown to be well approximated by a rigid rising 
spherical bubble of negligible velocity, leading to a rise velocity of (Davies and Taylor 1950) 



where g is the local gravitational acceleration and R is the radius of the spherical bubble. In 
the case of a two-dimensional bubble (eg, a rising cylinder), consideration of the potential 
flow over a cylinder rather than a sphere gives a numerical pref actor of 1/2 rather than 2/3. 

In the case of non-zero bubble density {e.g., non- unity Atwood number), the rising 
bubble will not 'feel' the full gravitational acceleration. Recently, a result has been derived 
(Goncharov 2002) for a related case, that of the Rayleigh- Taylor instability, which has been 
used in this context by lapichino et al. (2006): 



The combination of the terms 2At/(At + 1)5' can be thought of as an effective gravity, where 
2At/(At + 1) = 5p/pambicnt- Howevcr, while this correctly takes into account the density 
ratio, because of the different geometry being considered in the Rayleigh- Taylor instability. 




(7) 
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the numerical factor is inapplicable to the rising bubble case (that is, in the limit At ^ 1, 
this expression is off by a factor of (2/3)-\/27r ^ 1.67). 

A simpler expression one sometimes sees in the literature takes the form of the Davies 
and Taylor results and sets g Atg 



Vb. 



,DT+A = VAtgR . (9) 



This fails to capture the bubble motion in two respects — the incorrect At — ^ 1 limit, and an 
incorrect scaling with At. 

A further correction is due to the fact that our bubbles are expanding as they rise. 
Taking this into account {e.g., Ohl et al. 2003), together with the above correction for 
effective gravity on the bubble, we have 



/ 

d+1 

^6, rise " 



6 



2At {d-l)R\ {d-l)R 



-qR+ , 

At + 1^ \ 4 / 4 



(10) 



Here d = 3 for a spherical bubble and d = 2 for a 2d flame 'cylinder', and R is given by Eq. 6. 
Note that, all else being equal, expansion of the bubble acts to reduce the rising terminal 
velocity. 

As we see from the -R-dependence of Eq. 10, for larger and larger bubbles, the buoyancy 
increasingly outweighs any fluid drag on the bubble, and rise velocity grows. When the rise 
velocity exceeds the speed of the laminar flame, we would expect the shear on the sides and 
the vorticity generation behind the bubble to start to play a role in the bubble evolution. In 
particular, the bubble should begin to roll up. We note that this is also the point where the 
Rayleigh- Taylor instability will set in, as equating the bubble rise speed to the laminar speed 
yields the fire-polishing length (Timmes and Woosley 1992), although typically large-scale 
bubble distortions due to the bubble's own motion will occur first (Robinson et al. 2004). 

These bubble rise velocities are plotted a a function of the bubble sizes in Fig. 1 for the 
specific case oi a. d = 2 cylindrical bubble burning into a material of p = 4 x lO^g cm~^ in a 
gravitational field oi g = 10^'^ cm s~^, along with results of a multidimensional hydrodynam- 
ical simulation (described in §3) for comparison. We use Vb^nse as the fiducial laminar bubble 
rise velocity for the rest of this paper. 

In a turbulent medium, another source of drag exists beyond that created in the wake 
of the rising bubble — a 'turbulent eddy viscosity'. In the regime where this 'viscosity' domi- 
nates, the rise velocity coming from balancing the viscous drag for a sphere and the buoyant 
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force gives (Moore 1959) 

^^b,visc = ^ = TTT (11) 

t'turb -L^ fturb,6 

where the turbulent eddy viscosity on the scale of the bubble is assumed to be fturb = 
2i?Wturb,fe where Wturb,fe is the turbulent velocity on the bubble scale. We will use the same 
correction for effective gravity as above. 

Assuming Kolmogorov turbulence, for which there is some evidence in this context 
(Zingale et al. 2005b; Cabot and Cook 2006), the characteristic turbulent velocity on the 
bubble scale, (2/2), is 

Vtnrh,b = V[—] (12) 



L J 

where V is the turbulent intensity on the integral scale L. The corresponding velocity on 
the flame scale is 

//A //. 

VtnrhJ = V I j-j = ( ^ ) ^turb,fe. (13) 

Throughout this work we will assume characteristic values for V ~ 100 km on scales 
L ^ 200 km, consistent with, e.g. Hofiich and Stein (2002) and Kuhlen et al. (2006). 

In general, we expect the turbulent 'viscous' drag on the bubble to dominate for suffi- 
ciently small bubbles, and the drag produced by wake of the bubble itself to dominate for 
larger bubbles. The 'viscous' drag dominates if the viscous drag velocity prevents the bubble 
from rising as it would without the turbulent viscosity, f fe,visc < "^b.risc- Ignoring, for the time 
being, the effect of flame speed driven expansion on the rise velocity of the bubbles, this 
happens on scales 

/At + 1\^ (WlOO km 8-1)6 , , 

R < 2621.44 cm ^ . , ,n ,,,, , ^ , — 14 

V 2At J (^/lOio cm s-2)3(L/200 km)^ ^ ^ 

which, for typical values of At of 0.11, gives i? ^ 3.4 km. 

As long as the flame speed dominates the rise velocities, the evolution of the flame 
bubble will consist solely of the spherical expansion of the flame as smaller disturbances are 
burned out. As the bubble grows, the rise velocity also increases, leading to distortion of the 
sphericity of the bubble; a 'mushroom cap' forms caused by the vortical motions induced by 
the bubbles own rise. We refer to the radius at which the flame speed is no longer dominant 
as the distortion radius, R^. 

Beyond its effect on the rise velocity, the turbulent velocity field in the white dwarf on 
scales equal to or lesser than the bubble size can also directly distort the bubble. (In the 
extreme case that even the turbulent velocity on the scale of the flame thickness exceeds the 
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flame velocity, it can directly affect the burning, leading to the transition to the distributed 
burning regime (Niemeyer and Woosley 1997; Niemeyer and Kerstein 1997; Bell et al. 2004)). 

If the velocity of turbulent eddies of the size of the bubble is smaller than the flame 
speed, then the flame will burn through these, and little deformation will occur. As with 
buoyancy, however, as the bubble increases in size, the turbulent velocities on the scale of the 
bubble become larger and larger, and eventually they will begin to affect the dynamics of the 
bubble. This only happens on larger scales, however; because in the case we are considering 
the flame velocity is larger than the buoyant velocity, the turbulence can effect the slower 
buoyant rise speed before it can affect the faster-moving geometry of the flame bubble itself. 
Thus the distortion due to buoyant motions will occur first, and we need not consider direct 
turbulent distortion of the flame bubble here. 

Equating the viscous bubble terminal velocity in Eq. 11 with the laminar burning speed 
of the flame bubble in Eq. 6, we arrive at an expression for the maximum radius of a bubble 
before deformation sets in due to buoyant motions 

. 3.9 X 10' cm ( fWl°°k.ns-)3(V100Ums-n 

\^2At(l-At)y V (^/IQio cm s-2)3(L/200 km) J ^ ' 



oi Rk. 6.2 km for At = 0.11. We can do a little better than this by examining the behavior 
as a function of density. For the white dwarf model considered here, gravity is well fit (to 
within 5% between densities of 5 x 10*^ g cm~^ and 2 x 10^ g cm~'^) by 



g{p) = 1.198 X 10^° cm s~V(2-6269 - pg){p9 + 0.860752). (16) 

The properties of the flame from Timmes and Woosley (1992) are given in a fit provided in 
that paper, 

/ N 0.805 / ^/l2pN\ 0.889 

5, = 92 X 10^ cm s-M ^ ,) f^i-^) (17) 

however, we used this modified version of that flame speed which fits better at lower densities 
at the cost of some accuracy at higher densities: 



^■4(-(?)°"X^)°'''--(f)"l^^"'^'-- 

and we compute a similar fit for the Atwood number computed from the data tabulated in 
that paper. 

At = 0.114p90-3233x(12C)0-0479 (ig) 

where pg = p/lO^g cm~^ and X(^^C) is the mass fraction of carbon. 
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Substituting these results into Eq. 15, and expanding in a power series in log(p9) and 
log(X(^^C)) around a relevant density pg = 1.5 and X(^^C) = 0.3 where all of the fits are 
particularly good, and neglecting an additional very weak dependence on X(^^C), one finds 

R . 32^5 cm m''' f^V" ( ^ ( ^V". (20) 

If one instead considers the self-drag rather than the viscous drag, and follows the same 
procedure but with V given by Eq. 10 one finds 

Figure 3 shows this maximum radius where the 'viscous' rise speed becomes equal to the 
flame speed as a function of density, both from semi-analytically comparing the velocities 
as in Fig. 2, and from the expression given in Eq. 20. The plot also shows the transition 
from the rise velocity being dominated drag by turbulent eddies to that of the bubbles own 
motion at p ^ 2 x 10^ g cm~^. 

For concreteness, we consider these velocities given in Eqs. 6, 10, 11, 12, and 13 with the 
gravitational acceleration as a function of density taken from a pre-ignition Chandrasekhar 
white dwarf model produced with the Kepler code (Weaver et al. 1978) with a central density 
of 2.6 X 10^ g cm~^, a central temperature of 7 x 10^ K, and a carbon abundance by mass of 
0.3. Since the flame thickness varies considerably over the range of densities we examine, we 
measure the bubble radius in units of the laminar flame thermal thickness. The laminar flame 
properties — flame speed and Atwood number — are taken from Timmes and Woosley (1992), 
and logarithmically interpolated between tabulated points. Figure 2 shows the intersections 
of these different physical processes for bubbles with radii of 5, 500, 5 x 10^, and 5 x 10^ If. 

As we see, for small bubbles (R = 5 //), the laminar burning speed dominates over the 
turbulence and the buoyancy for all densities greater than ~ 8 x 10^ g cm~^. As the bubble 
gets bigger, the density below which laminar burning no longer dominates increases. For 
a bubble of radius 500 //, the crossover point is about 1.5 x 10^ g cm~^. For a 5 x 10^ If 
bubble, it is 3 X 10^ g cm"^, and for a 5 x 10^ // bubble, it is 8 x 10^ g cm~^. 

Clearly, at high densities, such as those of the ignition conditions, the bubbles can grow 
to enormous sizes (compared to the flame width) without feeling the effects of turbulence or 
buoyancy. As the panel for a 5 x 10^ // bubble shows, once the bubble is large enough, it is 
the buoyancy effects that will be felt first. This further suggests understanding the dynamics 
of flame bubbles is important to understanding the propagation of the very first flames in 
SNe la. 
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The previous discussion assumes a quasi-static flame bubble, considered at different 
sizes and positions/densities. The dynamic case, where the bubble expands and rises, is 
in principle much more complicated, but if we restrict ourselves to considering where the 
bubble is dominated by flame expansion, the picture greatly simplifles, as we see in Fig. 4. 
Here, we look at the evolution of individual bubbles with a variety of starting points and 
'evolve' the flame bubble in simple semi-analytic way — the size increases with laminar burn- 
ing speed (Eq. 6 and the position increases with buoyancy, assumed (as an upper limit) to 
be the minimum of the terminal rise velocities given above in Eqs. 11 and 10. We stop 
evolution when buoyancy speed exceeds flame speed, where we expect the bubble to start 
becoming distorted. As we see, because (by construction) the bubbles rise velocity is small, 
the 'evolution' consists simply of the bubble expanding to the critical size where rise velocity 
begins to compete with the flame velocity. Thus one can greatly simplify consideration of 
rising flame bubbles in high-density regions by ignoring hydrodynamical effects until they 
reach the size given by Eq. 20. A similar analysis was explored in Garcia-Senz and Woosley 
(1995) for igniting hot spots. 

Whether or not a flame burns through the center of the white dwarf can determine 
the importance of an off-center ignition {e.g., Plewa et al. 2004; Roepke et al. 2006), and 
can determine whether a central 'pool' of unburned fuel remains for possible later burning. 
Certainly if the undistorted bubble can expand to a size greater than its initial radial position 
from the center of the star, then it will both consume the fuel center of the star and proceed 
in more symmetric fashion. In Fig. 5, we consider the evolution of bubbles that initially begin 
with radii of 10 flame thicknesses and burn outwards. If the bubble begins sufficiently close 
to the center of the star — for the density structure of the progenitor considered here, within 
^ 23.5km — the flame will pass through the center. Woosley (2001) used similar arguments 
to show that the minimum separation of bubbles so that they do not merge as they float is 
about 5 km. Woosley (2005) used a simpler estimate of the terminal velocity of the bubble to 
estimate that the radius beyond which bubbles do not burn through the center is 20-40 km. 
This radius is sensitive to the central density of the white dwarf. 



3. NUMERICAL MODEL 

To test the ideas proposed above — and in particular the quantitative predictions of 
laminar expansion until buoyancy effects become important — we perform hydrodynamical 
simulations of burning, buoyant, flame bubbles. We use a low Mach hydrodynamics code 
(Bell et al. 2004) to follow the evolution. Here, the sound waves are flltered out of the 
system, allowing us to take much larger timesteps than a fully compressible code. A second- 
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order accurate approximate projection method is used to evolve the state. This code has 
been apphed to both terrestrial (Bell et al. 2005) and astrophysical flames (Bell et al. 2004; 
Zingale et al. 2005b). Since the laminar flames we are dealing with all have Mach numbers 
< 10~^, this is a very good approximation. 

We resolve the thermal thickness of the flame and track the dynamics of the bubble 
without the need for a flame model. Resolution studies with this code have shown us that 
we need about 5 zones in the thermal thickness of the flame to accurately capture the 
burning. A single ^^C C reaction rate is used, with the energy release corresponding 
to ^^Mg production. The rate is modifled to include screening (using the screening routine 
from the Kepler code, Weaver et al. 1978) — this leads to slightly higher flame speeds than 
reported in earlier works (see e.g. Bell et al. 2004). Finally, we use a general equation of 
state allowing for arbitrary degeneracy and degrees of relativity (Timmes and Swesty 2000), 
and the conductivities described in Timmes (2000). Table 1 lists the properties of planar 
laminar flames with this input physics. Resolving the flame restricts us to bubbles of small 
size (< 1000 //). However, by following bubbles in the flamelet regime, we expect that 
the results we flnd will scale to bubbles of large size, with similar Froude numbers. Our 
boundary conditions are outflow on the top and sides and slip-wall on the bottom. We note 
that these calculations are only used to test our ideas about the deformation radius of a 
flame bubble — we do not use these speciflc flame parameters in the semi-analytic estimates 
presented above. There, we use only the results presented in Timmes and Woosley (1992). 

We will follow the evolution of bubbles at different densities — and thus with different 
balances between flame speed and buoyancy — to get a feel for how the evolution unfolds and 
verify the analytic description sketched out above. The simulations presented here are all 
two-dimensional, using a Cartesian geometry. Previously, we looked at three-dimensional 
flame bubbles (Zingale et al. 2005a), but the range of length scales we need to resolve 
for a range of densities make a three-dimensional study infeasible. While axisymmetric 
calculations would capture the three-dimensional volume factors, the Cartesian geometry 
used here should give qualitatively the same behavior to the point of roll-up {e.g., Robinson 
et al. 2004). We do not consider externally-imposed turbulence in these simulations, as at 
these densities and scales the turbulence would swamp the effect we are trying to measure. 

Figure 6 shows the expected cross over points for where we expect buoyancy to dominate 
over the laminar burning for three fuel densities, 4.0, 3.0, and 2.35 x 10^ g cm~'^, assuming 
50% carbon / 50% oxygen. These plots were generated using the same method as Fig. 2, but 
considering flxed density and gravity and varying the flame bubble size. Flame properties 
(propagation speed and Atwood numbers) were interpolated from the values tabulated in 
Table 1. Turbulent velocities are shown for reference, but only the terminal rise velocity from 
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the bubbles own motions are shown. In these plots, as in Fig. 2, the laminar flame speed 
includes the Atwood number correction from Eq. 3 and the curvature correction from Eq. 4; 
because we are looking at larger flame thicknesses at these lower densities, the curvature 
correction is apparent in the variation of flame speed with bubble size. 

Figure 7 shows the evolution of a flame bubble with a fuel density of 4 x 10^ g cm~^. With 
our input physics, the laminar flame speed is 7.64 x 10^ cm s~^ and the thermal thickness is 
0.039 cm. The calculation was started by mapping a steady state laminar flame in a circular 
fashion onto the grid, with an initial radius of 6.4 If. Adaptive mesh reflnement was used 
to follow the flame. The flnest level of reflnement corresponds to a grid of 16384 x 32768 
zones (or 5 zones /If)- The domain is 128 x 256 cm in physical extent. As the flgure shows, 
initially, the flame simply burns outward, staying roughly circular. Only when it reaches 
a radius of ~ 500 // do we see signiflcant deformation. At this point, the bubble would 
continue to roll up, but we stop the calculation, since the bubble is starting to take up a 
considerable fraction of the domain width. 

The position of the bubble, its radius, rise velocity, and ash mass are shown in Figure 8. 
The top left panel shows the location of the topmost and bottommost extent of the bubble 
and the center of mass. The extrema are computed by laterally averaging the ash mass 
fraction and looking for the height where it exceeds 0.001. The center of mass is computed 
simply as 

E,,l/.A.^..^(^^Mg),, 

Here, the sum is over all the zones in the domain. Aij is the area of the zone {Ax Ay) , pi j is 
the mass density of that zone, X{'^^Mg)ij is the mass fraction of ash in the zone, and yj is 
the vertical coordinate of the zone. The top right panel shows the radius of the bubble. This 
is taken to be half the distance between the left and right extrema of the ash, computed 
as described above. This is more properly a width of the bubble. This will be the true 
radius for the initial evolution of the bubble, when it is circular. As the bubble deforms, 
this is actually a measure of the cross-sectional radius of the bubble. The bottom left panel 
shows the velocity of the center of mass at timelevel n, v"'. This is computed via centered 
differencing of the center of mass position at the points we output plot flies. Finally, the ash 
mass, shown in the bottom right panel, is computed as 

M = J2p^,JA,Xe'Mg),, (23) 

id 

Since these are two-dimensional simulations, this has units of g cm~^. 

As we see, for the flrst 10~^ s, the bubble is expanding, and the downward flame prop- 
agation is faster than the rise velocity of the bubble. Then the bubble rise velocity becomes 
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equal to the laminar flame speed, and the ash region no longer reaches lower heights. For 
most of the evolution of the bubble, the radius expands according to Eq. 3, as expected 
when the buoyancy effects are negligible. The velocity of the bubble steadily rises through 
the course of the simulation, staying very close to the predicted terminal velocity (Eq. [10]) 
for the initial evolution. Since the bubble starts from rest and is continually accelerating, 
the velocity curve should not lie directly on top of the terminal velocity curve. Only at very 
late times, when the bubble is no longer spherical, does the rise velocity from the simulation 
overtaking our prediction. The mass plot shows that the total mass of the bubble increases 
by 4500 X over the simulation, showing the dominance of the burning at this density. We 
stop the simulation once the bubble grows to more than half the width of the domain to 
prevent any influence from the boundaries. 

At 3 X 10^ g cm~^, the evolution proceeds in much the same fashion (Figure 9). Because 
the flame speed is lower {Si = 4.34 x 10^ cm s~^), we expect that the bubble will begin to 
roll up at a smaller radius. This simulation was run in a domain 192 x 256 cm large, with an 
effective grid of 12288 x 16384 zones (or 5 zones / //). The steady state flame was mapped 
onto the grid with an initial radius of 4.1 As the figure shows, the bubble burns radially 
outward until it grows to a radius of ~ 100 //. At this point, it begins to deform. The 
integral quantities (Figure 10 again show the flame initially propagating downward until the 
rise velocity exceeds the laminar flame speed. We also again see that our prediction of the 
bubble velocity again agrees with the simulation to the point where the bubble begins to roll 
up. 

Figure 11 shows the evolution of our lowest density bubble (fuel density of 2.35 x 
10'' g cm~'^). The laminar flame speed is roughly 3x lower {Si = 2.58 x 10^ g cm~^) and the 
flame thickness is If = 0.18 cm. The domain size is 160 x 320 cm, with an effective grid of 
4096 X 8192 zones. The flame solution was mapped onto the grid with an initial radius of 
8.8//. As discussed above, we expect this bubble to roll up at a much smaller radius. As the 
figure shows, already at a radius of 30//, the bubble is deforming. 

The position, radius, velocity, and mass for the 2.35 x 10^ g cm~^ bubble are shown 
in Figure 12. As we see, the bubble quickly inflates in the very earliest moments, as the 
planar flame used to initialize the flame settles into the circular configuration. For the first 
5 — 6 X 10~^ seconds of evolution, the radius evolution parallels the theoretical prediction. It 
then gradual begins to diverge as the bubble begins to roll up. The buoyant velocity of the 
bubble quickly exceeds the prediction (Eq [10]), but levels off and stays roughly 1.25 — 1.5x 
higher for the remainder of the evolution. The mass increases by almost 120 x through the 
course of evolution. 

The low density calculation was run at a resolution corresponding to 4.6 zones in the 
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thermal flame thickness. To test the sensitivity of our results, we performed a convergence 
study. Figure 13 shows the bubble velocity and mass as a function of time for four different 
resolutions, corresponding to Zj/Ax of 1.15, 2.3, 4.6, and 9.2. As we see, the three highest 
resolution curves track each other very well, indicating that our simulations have converged. 
Only the very lowest resolution run shows a significant divergence, producing about 20% less 
ash than the converged solution. This is the result of under resolving the flame. 

Figure 14 shows the carbon destruction rate {—dX(^'^C)/dt) for the 2.35 x 10^ g cm~^ 
bubble at 10~^ s for the high-resolution run. A planar laminar flame has a peak carbon 
destruction rate of —1.3 x 10^ s~^. We see that the shear at the top of the bubble suppresses 
the burning rate slightly, while the cusp near the bottom of the bubble has a significantly 
higher burning rate. The flow field is represented by the vectors superposed on the rate. 
We see than below the bubble, the flow is carrying fuel toward the flame while above it, 
the flow is carrying fuel away from the flame. This flow pattern is similar to that seen in 
Zingale et al. (2005a), where a lower density flame bubble was evolved in three-dimensions. 
There, the flow managed to punch a hole through the top of the bubble, transforming it into 
a torus. 

To within constants of order unity, the cross over point predicted in the analysis in § 2 
seems to agree with the numerical results. In our two higher density bubbles, we had an 
extended period where the bubble stayed roughly spherical, and the rise velocity showed 
excellent agreement with our theoretical prediction, Eq. 10. This suggests that this picture 
can safely be extended to the ignition conditions, to allow us to make predictions about the 
evolution of the very first flame bubbles. 

4. DISCUSSION 

We presented a model for the early evolution of flame bubbles. Several functional forms 
for the terminal velocity of flame bubbles were explored. Our analysis was supported by 
resolved numerical simulations of flame bubbles at a range of densities. We showed that for 
initial flame bubbles at densities of pg > 2 (eg, the inner ^ 250km of the white dwarf) the 
laminar flame speed dominates other motions at least until the flame bubble approaches on 
order a half kilometer in size. This distance from the center is comparable to that explored in 
a recent study of off-center ignition (Roepke et al. 2006). A bubble at the very center of the 
star (pg = 2.6) will deform once it reaches a radius of 1 km. This suggests an (approachable) 
minimum resolution to capture all of the early bubble dynamics in large-scale simulations 
is on the order of a tenth of a kilometer. For comparison, the current state of the art 
explosion models (e.g. Ropke and Hillebrandt 2005) seed initial flame bubbles with radii 
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of a few kilometers with a resolution of 1 km, so they are just able to resolve the relevant 
lengthscales for a bubble ignited precisely at the center (by putting only 2 zones across the 
entire bubble). 

These results further suggest a minimum size to which turbulent fragmentation can 
shred a burning region. For instance, in lapichino et al. (2006), where a non-burning bubble 
is simulated, the turbulent motions dominate and it is suggested that this will continue to the 
smallest scales on which a flame can be ignited. However, the results presented here imply 
that had flame physics been included in those simulations the cross-over point between the 
turbulent speeds and flame speeds would have been well-resolved, and fragmentation would 
have ended on the scale of hundreds of meters, rather than tens of centimeters. 

The existence of a (high) minimum scale to fragmentation places a limit on the burning 
rate of the initial rising burning region. If a burning region of volume V is continuously 
fragmented into near-spheres of radius Rf, then the total surface area available for burning 
is 3V/Rf and the burning rate is 

dV /3V~ 



Using our simple fits for Rf{p) and Si{p), neglecting the At number effect as small < 25% 
at high densities, and assuming that we are dealing with large enough regions that flame 
curvature effects are small, we have a burning law 

H . l.a4 X 10' =ec- - ( y \ (^) (^) ^° ,25) 
V ^ VlOOkms-V V200 kmy V 0.3 7 ^ ' 

Note that this is a different picture from the usual 'turbulent flamelet' sort of burning 
model, based on simulations such as Khokhlov (1995). In that paper, fragmentation was 
suppressed both by the geometry (a planar flame front, externally imposed by symmetry, 
so that on the largest modes fragmentation was suppressed) and the extremely thick model 
flame (so that only a few flame thicknesses fit in the box, and only the largest scales if even 
those could fragment). 

In the turbulent flamelet picture, it is explicitly assumed there is a large separation of 
scales between that of the buoyant plume and the flame thickness, so that the flame can be 
greatly wrinkled without altering the overall geometry of the burning. This certainly seems 
to be a reasonable picture at large- to intermediate- densities (p ^ 5 x 10^ g cm"'^) where 
the typical fragment size would be tens or hundreds of thousands of flame thicknesses. 

However, in non-planar, 'rising bubble' geometries it is seen that if a bubble can frag- 
ment, it does, on roughly the time it takes to rise one bubble height (Robinson et al. 2004; 
lapichino et al. 2006). At high densities this would happen very slowly, as bubbles are large 
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and rise velocities are small. However, considering Fig. 3, we see that for lower densities 
(p < 2 X lO^g cm~^), the typical fragment size would be hundreds of flame thicknesses or 
less, in which case fragmentation would quickly outstrip wrinkling as a method of produc- 
ing burning surface area, producing a rapid increase in burning at an interesting density 
for reproducing observables in Type la supernovae. This will all happen on scales much 
smaller than those resolvable by multidimensional simulations. Verifying this new burning 
mechanism will require three-dimensional simulations to investigate fragmentation. 

Another refinement to our model would be to include some stochastic advection, which 
changes not only the flame position, but more importantly the burning conditions — flame 
speed and local gravity. In particular, this would affect the radius at which ignited bubbles 
are guaranteed to burn through the center. Numerical simulations of the convective phase 
of SNe la have shown that large scale motions can exist. These can carry the flame bubbles 
far away from their initial location, affecting the results presented here. 
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Table 1. Properties of planar simplified ^^C + ^^C ^^Mg flames into given fuel 

conditions. 





X(i2C) 


Si 




At^ 




I cm^'^ 


cm 


cm 




2.35 


0.5 


2.58 X 10^ 


0.18 


0.241 


3.0 


0.5 


4.34 X 10"^ 


0.078 


0.224 


4.0 


0.5 


7.63 X 10"^ 


0.039 


0.205 



^Density of fuel in units of lO^g cm ^. 

'^Flame thermal thickness, If = (T^sh — 
Tfuci)/ max{VT}. 

'^At = (pfuel - Pash)/ (Pfucl + Pash)- 
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Fig. 1. — Different parameterizations of predicted rise velocity of a buoyant rigid spherical 
bubble (lines), for the specific case of a cylindrical flame burning into material of density 
p = 4 X lO^g cm""^. Plus symbols represent simulation results described in §3. Lines are 
plotted for, top to bottom: dashed line - Davies- Taylor result, w;,,dT) Eq- 7; dotted line, same 
functional form, but with simple Atwood-number correction, fb,DT+A; Eq- 9; thin solid line, 
rise velocity fb^rise from Eq. 10 but not taking into account flame expansion; solid line, rise 
velocity Vh^rise from Eq. 10 taking into account flame expansion; dot-dashed line, Goncharov 
result t'ft.RT from Eq. 8 for Rayleigh- Taylor instability. All of these plots are for terminal 
velocities, so provide only an upper envelope for bubble speed. Horizontal line represents 
laminar planar flame speed, and vertical line indicates where buoyancy becomes dominant 
(here defined where the bottom of the bubble begins moving upwards) and so bubble begins 
deforming and these velocities no longer strictly apply. 
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Fig. 2. — Velocities affecting the propagation of a flame bubble of given initial size at various 
densities in our white dwarf model. Plotted are the flame expansion velocity R (solid line), 
fiducial turbulent velocities (dashed lines) assuming large-scale convective motions on 200km 
scales at ~ 100km s~^; turbulent velocities are calculated on scales of the bubble (fturb.fe, 
upper dashed line) and flame thickness (t'turb,/, lower dashed line). Also plotted are terminal 
rise velocities of the flame bubble assuming the bubble rise speed is dominated by drag 
induced by the bubbles own motion (t^b^nso, + symbols) or determined by turbulent drag 
(^^'6,visc, o symbols). These quantities are plotted, left to right and top to bottom, for flame 
bubbles of size 5, 5 x 10^, 5 x 10^, and 5 x 10^ flame thicknesses at conditions corresponding 
to the given density. Buoyant and turbulent velocities are comparable because the origin of 
the large-scale turbulence is buoyant convection. 
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Fig. 3. — Plot in (p, r) space showing the size of a bubble r (on left, in local flame thicknesses; 
on right, in cm) where laminar flame speed dominates the other speeds calculated here, using 
same model as in Fig 2. Plus symbols indicate calculated values, and the solid line indicates 
fit values: o„ left R = Sll,pl<"\ on right R = 12,67cmpi--. On the right is also plotted 
(the top solid line) the analytic approximation for the distortion radius for a turbulent rise 
velocity as given in Eq. 20 (dotted) and that for the laminar rise velocity (dashed). 
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Fig. 4. — Semi-analytic evolution of flame bubble sizes and positions, shown as trajectories in 
{R, r) and (p, r) space, with flame bubbles of varying sizes and initial positions. The top panel 
shows the evolution of bubble radius r and position from the center of the white dwarf R; the 
bottom panel shows position in terms of density p. Trajectories are evolved from the initial 
conditions (o) until the flame velocity no longer dominates and evolution becomes turbulent 
and fragmentation of the bubble becomes possible (marked by 'X'). Varying thickness of 
the lines denote different initial conditions, with lines of the same thickness in each panel 
indicating evolution of the same bubble in the different plots. Because by definition we 
are only considering times where while buoyant rise velocities are small compared to flame 
velocities, the position and ambient density change little during the course of the evolution, 
although at very high densities where the flame speed is quite high more rising can occur. 
The bubbles are initially set to be 10 and 5000 flame thicknesses in radius; as this figure 
shows, the only difference in evolution in these cases is where along the trajectory the flame 
bubble begins. On the bottom panel, the dotted line indicates the simple fit for the critical 
bubble size as a function of density, as plotted in Fig 3, and the dashed line indicates the 
analytic approximation given in Eq. 20. 
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Fig. 5. — Semi-analytic evolution of flame bubble sizes and positions, as in Fig. 4, but 
focused on bubbles which start very near the center - plotted are trajectories of bubbles 
with initial positions, left to right, of 21, 22, 23, 24, 25, 30, and 50 km - and the end of 
the trajectories. Those bubbles whose final size before distortion by turbulence or buoyancy 
exceed their radial positions certainly and immediately burn through the center of the star, 
and thus do not leave a central region of fuel behind. The first bubble whose final size does 
not exceed their radial position sets the maximum size of a bubble which does not burn 
through the center, because of the downward and monotonic relationship between density 
and final size as shown in Fig. 3. For the white dwarf progenitor considered here, bubbles 
beginning at positions closer than ^ 23.5km to the center of the star will burn through the 
center, while those at a larger distance, while burning to a radius of ~ 30km before other 
effects take hold, will not. 
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Fig. 6. — Plot of predicted velocities for a flame bubble burning a 50/50 mixture by mass of 
Carbon and Oxygen, using the simplified ^^C + ^^C — > ^^Mg burning network, for comparison 
with numerical results. Shown are quantities at density p = 2.35 x lO^g cm~^, (top) with a 
flame thickness of 0.18 cm, p = 3 x 10''g cm~^, (middle) with a flame thickness of 0.078 cm, 
and p = 4 X lO'^g cm~^ (bottom) with flame thickness of 0.039 cm. Gravity is fixed at 
g = lO^'^cm s~^. The solid line shows laminar flame velocity, plus signs indicate terminal 
bubble rise velocity, and dashed lines indicate turbulent velocities on bubble (top) and flame 
(bottom) scales assuming motions of 100km s~^ on an integral scale of 200km. Terminal 
bubble speed is plotted assuming a quiescent medium, e.g. drag is due to the bubbles 
motions only. 
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Fig. 7. — Evolution of a 4 x 10'' g cm~^ flame bubble. The bubble burns out in a circular 
fashion until it reaches about 400 thermal thicknesses in radius. At this point, it begins to 
deform. Only the lower half the computational domain is shown here. 
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Fig. 8. — Bubble position, size, rise velocity, and mass as a function of time for the 4 x 
10^ g cm~^ flame bubble. In the top left panel, positions are shown for top and bottom of 
the bubble (solid lines) and center of mass of the bubble (dashed line). In the top right 
panel, radius is plotted, in both centimeters and flame thicknesses, with the solid line being 
the measured size and the dashed line given by Eq. 3. In the bottom left, the bubble rise 
velocity is plotted as measured (solid line) and as given by Eq. 10 (dashed). On the bottom 
right panel is shown the burned mass inside the bubble. 
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Fig. 9. — Evolution of a 3 x 10'' g cm~^ flame bubble. The bubble burns out in a circular 
fashion until it reaches about 100 thermal thicknesses in radius. At this point, it begins to 
deform. Only the lower portion of the computational domain is shown here. 
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Fig. 10. — Bubble position, size, rise velocity, and mass as a function of time for the 3 x 
10^ g cm~^ flame bubble. In the top left panel, positions are shown for top and bottom of 
the bubble (solid lines) and center of mass of the bubble (dashed line). In the top right 
panel, radius is plotted, in both centimeters and flame thicknesses, with the solid line being 
the measured size and the dashed line given by Eq. 3. In the bottom left, the bubble rise 
velocity is plotted as measured (solid line) and as given by Eq. 10 (dashed). On the bottom 
right panel is shown the burned mass inside the bubble. 
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Fig. 11. — Evolution of a 2.35 x 10^ g cm~^ flame bubble. In contrast to Figure 7, the 
lower density bubble begins to deform at a much smaller radius. Only a portion of the 
computational domain is shown here. 
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Fig. 12. — Bubble position, size, rise velocity, and mass as a function of time for the 2.35 x 
10^ g cm~^ flame bubble. In the top left panel, positions are shown for top and bottom of 
the bubble (solid lines) and center of mass of the bubble (dashed line). In the top right 
panel, radius is plotted, in both centimeters and flame thicknesses, with the solid line being 
the measured size and the dashed line given by Eq. 3. In the bottom left, the bubble rise 
velocity is plotted as measured (solid line) and as given by Eq. 10 (dashed). On the bottom 
right panel is shown the burned mass inside the bubble. 
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Fig. 13. — Bubble radius and rise velocity as a function of time for the 2.35 x 10^ g cm"'^ 
bubble, with simulations run at four different resolution. The solid line represents results 
using the resolution used in the present analysis (4.6 zones / //), the dotted line is a higher 
resolution run (9.2 zones / //), the dashed line is a low resolution run (2.3 zones / //), and 
the dot-dash line is at a lowest resolution (1.15 zones / //). The fiducial- and high-resolution 
curves lie nearly on top of one another, demonstrating that our simulation is converged. 
Significant divergence of the very low resolution run in the mass plot is seen — evidence that 
at this lowest resolution, we are incompletely resolving the burning, although the large-scale 
rising dynamics is robust even at this very low resolution. 
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Fig. 14. — Base-ten logarithm of the carbon destruction rate {—dX(^'^C)/dt) for the 2.35 x 
10^ g cm~^ bubble at 10~^ s, after the deformation has begun. For reference, a planar 
laminar flame has a peak destruction rate of 1.3 x 10^ s~^, or \ogiQ{—dX{^^C)/dt) = 5.1. 
Notice that the most vigorous burning occurs on the bottom of the bubble, with an intensity 
approximately twice that of an unstrained laminar flame. 



